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The major accomplishments during this period are summarized in the 
enclosed abstract which was submitted for the Aerospace Sciences meeting. 
Other approaches to turbulence modeling are being pursued. 
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Abstract 

A one equation turbulence model based on the turbulent kinetic energy 
equation is presented. The model is motivated by the success of the Johnson- 
King model and incorporates a number of features uncovered by Simpson’s 
experiments on separated flows. Based on the results obtained , the model 
duplicates the sucesses of algebraic models in attached flow regions and out 
preforms the two equation models in detached flow regions. 


Introduction 


*«*#••• 


This work is the first phase of an on going research effort designed to develop 
turbulence models for separated turbulent flows over transonic airfoils. It 
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is generally accepted that algebraic turbulence models are not suited to 
describe non-equilibrium turbulent flows. This was clearly brought out by 
the extensive work of Simpson and his associates (see reference 1 for a 
summary ). The next logical step is to use a higher order closure model. 
Unfortunately, the work of Patel, Rodi and Scheurer [2] showed that most of 
the fc — e models and related two equation models are unable to predict some 
simple flows. Moreover, the work of Coakley [3] showed that predictions 
of'two equation models are not much better than algebraic models in the 
presence of separation. 

The recent model of Johnson and King [4], which is based on a simplified 
version of the tubulent kinetic energy equation, preformed well in regions 
where non-equilibrium turbulent effects are important. However, it did 
not always preform well in regions where equilibrium turbulent flow exists. 
This negative result does not diminish the significance of their contribution. 
Rather, it suggests that more physics exists in the turbulent kinetic energy 
equation that is yet to be explored and exploited. 

Because of the above, an effort was undertaken to develop a one equa- 
tion turbulence model based on the turbulent kinetic energy equation for 
both attached and detached flows. Simpson made two relevant conclu- 
sions. First, scaling within and outside separated flow regions is not the 
same. Second, the eddy viscosity assumption is incorrect in the backflow 
region of a separating boundary layer. Thus, the model developed retains 
the accepted form of the turbulent kinetic energy equation but adjusts the 
•tubulent shear stress and the scaling depending on the local flow conditions. 

Most of the previous work using one equation models employed the 
same length scale for turbulent stresses and turbulent dissipation. This led 
to numerical difficulties in the near wall region with the result that most 
researchers lost interest in such a model. The difficulties can be avoided 
by using two different length scales in the near wall region, one for the 
’turbulent stresses and one for turbulent dissipation. Specification of the 
length scales requires close attention to scaling near the wall. 

The ,;N avier- Stokes code employed is based on that used by Swanson 
and Turkel [5]. The turbulent kinetic energy equation is coupled to the 
conservation equations and a four^stage Runge-Kutta time stepping scheme 
is employed. The Reynolds averaged N avier- Stokes equations that use the 
Johnson- King model or higher order closure models cannot be started from 



a uniform flow condition. The usual procedure is to use an initial solution 
obtained from a flow solver that incorporates an algebraic turbulence model. 
This innefflcient starting proceedure is not necessary in the present work. 


Approach 


The equations governing turbulent flow past airfoils are the Reynolds aver- 
aged conservation equations of mass, momentum and energy and the tur- 
bulent kinetic energy equation. When Favre’s mass averaging and indicial 
notation are used the conservation equations are ; 


P,t + {puj) ,j = 0 . 

(pui) , t 4- [pUjUi +SijP - ( fij - pu"u")] ,j = 0. 

(pE) ,t + [*i (pE + P) + qj -(- pu'jh" - u i (t^ - pu?«?)] ,j = 0 


where 


r_.2-.-i 

Tij — ft u iij "hit j,i g 


qj — AT V - 


P = 7-1) 



1 . _• 

2 UiUi . 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


where p is the density, U{ is the mean velocity in the direction of X{, E and 
h are the total energy and enthalpy per unit mass, P is the mean pressure, 
ft , A and T are respectively, the molecular viscosity and conductivity and 
temperature. In attached flow regions, the turbulent stress and heat flux 
are determined from the eddy viscosity (Bousinessq) approximation. 
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Pr t is the turbulent Prandtl number. fi t is the eddy viscosity and k is the 
mean turbulent kinetic energy per unit mass. 


1 


pk = -pu'yi 


(9) 


For detached flow the above, relations are replaced by correlations obtained 
from Simpson’s experimental data. 

The modeled turbulent kinetic energy equation is; 

(pfc) , t + pUjk - v ( pic ) ,j -c k — pu'yj ( pie ) 


„ = -pu'y-Uijj ~pt (io) 

where e is the turbulent kinetic energy dissipation rate. 

( 11 ) 




e = 




c k is a constant (0.1). i t is the dissipation length scale and u is the molecular 
kinematic viscosity. In a one equation turbulence model the eddy kinematic 
viscosity is given by 

u t = ( 12 ) 

where is a constant and is the turbulent length scale. 

Simpson emphasized the importance of normal turbulent stresses in 
separating flows regions. Because of this, the governing equations retain all 
normal stresses. 


The Length Scales 

To complete the formulation of the model, the length scales must be de- 
fined. The expressions are summarized here and a complete derivation 
will be given in the paper. All variables are non-dimensionalized with the 
chord being the characteristic length scale and {Pool Poo)* the characteristic 
velocity. Letting 

L. = C,y [l - exp 

Lu = C,y [l - exp (^)j 

n s') 
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then 


£ e = min(L t ,y mas ) 

t-H — min (1/^) Umax) 

where y^ a x is the value of y at which F(y) is maximum, ui is the vorticity 
and R e is the Reynolds number. The constants in the definition are 
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= 0.09 

K 

= 0.41 
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An 

= 76. 

Af 

= 47.47 


Results and Discussion 

The results presented here are for a flate plate at zero angle of attack and 
for am NACA 0012 airfoil at a variety of Mach and Reynolds numbers and 
angles of attack. Fig. 1 shows flat plate results for a free stream Mach 
number, M of 0.5, Re = 10® and an amgle of attack, a , of zero. As is 
seen in the figure, good agreement is indicated for the skin friction, Fig. la, 
velocity distribution in the near wall region, Fig. lb, amd the one seventh 
power law, Fig. lc. 

The next set of calculations are for an NACA 0012 aurfoil. All calcula- 
tions employed a 160 X 79 C-grid. The normal spacing of the first point 
off the wall was 2.XT0 -6 . The first case considered is that for M w = 0.5, 
Re = 2.89.XT0 6 and a = 0.. Fig. 2(a) compares the pressure distribution 
with the measurements of Ref. 6. Fig. 2(b) compares the skin friction 
coefficient with a boundary layer calculation given in Ref. 7. As is seen 
in the figure, good agreement is indicated. The difference in the skin fric- 
tion coefficients at the leading amd trailing edges aure a result of much finer 
resolution in the boundary layer calculation. Another comparison with the 
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experiment of Ref. 6 is given for M = 0.756 , a = 0. and = 4. 01 AT 0 6 . 
Again excellent agreement with experiment is obtained. 

All of the above cases are for a = 0. . Fig. 4 compares predictions of 
this theory with the experiment of Ref. 8 for a = 1.49° and R^ = 9.A10 6 .. 
As is seen in the figure, both upper and lower surface presure distributions 
are well predicted by the present one equation model. 

Simpson noted that in the backflow region of a separating boundary 
layer, the eddy viscosity assumption is incorrect. Because of this, Simpson’s 
measurements are used to develop correlations for the turbulent stresses in 
the backflow region. Unfortunately, there is scatter in the data . Bacause 
of this, the results presented in Fig. 5 ( M ’* = 0.799 , R* = 9.AT0 6 , 
a = 2.26° ) should be considered preliminary. 

Details of the separated models together with comparisons for the RAE 
2822 and the Cast 10 airfoils will be presented in the paper. Based on the 
results obtained to date, it is concluded that one equation models, properly 
formulated, are capable of out preforming two equation models at reduced 
computational cost. 
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Figure la. Turbulent flat plate flow ( M * 0.5, Re = 10 ) 

comparison of skin friction distributions. 



Figure lb. Turbulent flat plate flow ( M = 0,5, Re = 10** ): 
comparison of velocity profiles. 
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Figure 2b. Turbulent flow over an NACA 0012 airfoil ( M = 0.$ 
Re = 2.89xl0 6 ): skin friction distributions. 
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Figure 3. Turbulent flow over an NACA 0012 airfoil ( M = 0.756, 
di = 0 * , Re = 4.01xl0 6 ): pressure distributions. 
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Figure 5b. Turbulent flow over an NACA 00] 2 airfoil ( M = 0.7??, 
ot = 2.26* , Re = 9.x]0 6 ): pressure distributions 
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